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- - - Abstract 

I We present predictions for /fj-meson form factors obtained from the analysis of QCD sum 

rules in next-to-leading order of perturbation theory. The radiative corrections turn out to be 
sizeable and should be taken into account in rigorous theoretical analysis. 
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1 Introduction 



The method of QCD sum rules PP is designed to estimate low-energy characteristics of hadrons, such 
as masses, decay constants and form factors. Within this framework we analyze the correlation func- 
OO . tion of currents in deep euclidean region with the help of operator product expansion, which allows 
^ ', us to take into account both perturbative and nonperturbative contributions. The presence of latter 
could be traced to the non-vanishing values of vacuum QCD condensates. Physical quantities, we 
are interested in, are determined by matching this correlator to its phenomenological representation. 
O ' In this work we performed an analysis of three-point sum rules for p-meson form factors at 
^ intermediate momentum transfer. Basically, it is an extension of already available LO analysis pj to 
Ph! include radiative corrections. To compute radiative corrections we used the technic already developed 
PhI and tested in the analysis of pion electromagnetic form factor within NLO QCD sum rule setup both 
^ I with pseudoscalar and axial- vector pion interpolating currents 01^. 

The paper is organized as follows. In section 2 we describe our framework and give explicit 
expressions for next-to-leading order corrections to double spectral density. Section 3 contains our 
numerical analysis and expressions for the contributions of gluon and quark condensates. Finally, in 



■ section 4 we draw our conclusions. 



2 The method 

To determine p-meson electromagnetic form factors we will use the method of three-point QCD sum 
rules. Within this framework p-meson is described as a result of an action of vector interpolating 
current on vacuum state. We define the vacuum to p-meson transition matrix element of vector 
current as 

2 

in 

{0\j^\p+,e) = ^e„ (1) 
9p 

where nip is the p-meson mass, gp is the p — 7 coupling constant {gp/4:7i = 1.27) and stands 
for p-meson polarization vector. Next, assuming parity and time- reversal invariance, the general 
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expression for p-meson electromagnetic vertex could be written in terms of three form factors: 



+ [9^.aqp - 9f.pqa\F2{Q^) + \pMp + p )^F^{Q^) \ , (2) 

nip ) 

where j^' = euUjfj,u + edd'jfj^d, the momenta of initial and final state p-mesons were denoted by p, p' 
and = —q^{q = P — p') is the square of momentum transfer. Fi, F2 and F3 are electric, magnetic 
and quadrupole form factors. At zero momentum transfer these form factors are expressed through 
the usual static quantities of p-meson charge, magnetic moment (p) and quadrupole moment {D): 

Fi(0) = 1, (3) 
F2(0) = (4) 

m) = ^{i^-i)-^D. (5) 

Within approach of QCD sum rules the theoretical estimates of p-meson form factors follow from 
the analysis of the following three-point correlation function: 

ii,o.0 = t'J d'xd'ye^^^'--^-y\o\T{j;ix),j;\o),Uy)m, (6) 




Figure 1: LO diagram 



The scalar amplitudes Ilj in front of different Lorentz structures are the functions of kinematical 
invariants, i.e. Ilj = nj(p^,p'^, g^). In the region of large euclidean momenta j9^,p'^,g^ < this 
correlation function could be studied with the use of ordinary perturbative QCD. The calculation 
of QCD expression for three-point correlator is done through the use of operator product expansion 
(OPE) for the T-ordered product of currents. As a result of OPE one obtains besides leading 
perturbative contribution also power corrections, given by vacuum QCD condensates. We will return 
to the discussion of power correction to QCD sum rules later after the definition of Borel transform 
for our correlation function. Now let us discuss the perturbative contribution. The calculation 
of perturbative contribution could be conveniently performed with the use of double dispersion 
representation in variables Si = p^ and S2 = p'^ at < 0: 

n;::;;*(p',p", = ^ j J^fz^J^^rJ^^^'^"^'^ + subtractions (7) 
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The integration region in ((7|) is determined by condition ^ 

-1< <i 



(8) 



and 



A(si, X2, 0:3) = (xi + a;2 - Xa)^ - 4xiX2. 



(9) 



The double spectral density p^^i^isi, S2, Q^) is searched in the form of expansion in strong coupling 
constant: 



PlZ'is,, .2, Q') = P^^s,, .2, Q') + (^) s„ Q') + 



(10) 






Figure 2: NLO diagrams 



At leading order in coupling constant we have only one diagram depicted in Fig. 1, contributing 
to three-point correlation function. At next to leading order we have 6 diagrams shown in Fig. 
2. The calculation of corresponding double spectral density was performed with the standard use 
of Cutkosky rules. In the kinematic region < 0, we are interested in, there is no problem in 
applying Cutkosky rules for determination of p^aisisi, S2,Q'^) and integration limits in si and S2- 
The non-Landau type singularities, not accounted for by Cutkosky prescription, do not show up 
here. The calculation could be considerably simplified with the use of Lorentz decomposition of 
double spectral density based on a fact, that our spectral density is subject to three transversality 
conditions: p^a^Qf, = Pf^afiPa = Pp^apP'p = 

= A, [{Q' + - (x + y)p^] [{y - x)p^, + {Q' + ] \{Q' + y)pt + {Q' - y)p^] 

-\a,[{q' + x)pi - (x + yK][2(pf - K + {Q' + y)9'^] 

-\a,[{x - - (Q' + [2(P2 - + {y- Q')9'% (11) 

^In our case this inequality is satisfied identically. 
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where x = S1 + S2, y = S1 — S2 and pi = p,P2 = p' ■ The four independent structures Ai (we suppressed 
the dependence on kinematical invariants) are given by a solution of system of hnear equations: 

h = p^./3Pi'p^P? = y(A;Ai-A2-A3-A4) (12) 

h = p^af3Pt9''^ = jix + Q'')(^kA,-3A2-As-A,^ (13) 

h = p^a/3P2 9'''='liy + Q'){kAi-A2~3A,-A,^ (14) 

h = p^a/3pfg'"' = ~{y-Q')[kA,-A2-A3-3A,^, (15) 

where k = A(si, S2, —Q"^)- Solving this system it is easy to find explicit expressions for Ai in terms 
of li (functional dependence on kinematical invariants is assumed): 

20 2 2 2 

- F^-(g^^- 

At Born level and expressions for are easy to find and they are given by (S3 = Q"^) |2]: 

^1 - 2A;V2 ' 

(0) _ _3£i£2 , . 

^2 - ^1/2 ' \'^^) 

= (22) 



(0) _ _3£iS3 

-'4 - ^1/2 ' y^'^i 

The calculation of NLO radiative corrections to double spectral density is in principle straightfor- 
ward. One just needs to consider all possible double cuts of diagrams, shown in Fig. 2. However, 
the presence of coUinear and soft infrared divergences calls for appropriate regularization of arising 
divergences at intermediate steps of calculation and makes the whole analytical calculation quite 
involved. We will present the details of NLO calculation in one of our future publications. Here we 
give only final results: 

^1/2 jW ^ __g3 _^ s^Si + S2S1 - S2 + (Si + 32)^3 - S^Ss - S2S3 + S1S2S3 

-161og^(i;i) - 161og(f3) log(fi) - 161og(t;4) log(t;i) 

+2 log(t;i) - 4 \og^{v^) - 4 W^v^) - 2 log(t;2) - 2 log(t;3) - 8 log(t»3) ^og^v^) 

(I) - (I) - (t) - (I) ^ (1)1 • <^*> 

A;l/2ja) ^ _2s2 _ 2^2 ^ 2S^ - 8S1S2 + S1S2 
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-32 - 32 log(^;3) log(^;i) - 32 log(^;4) log(^;i) + 4 log(vi) 

-81og^(^;3) - 81og^(T;4) - 4:log{v2) - 4:log{vs) - 161og(T;3)(T;4) 



X2\ 



XiJ 



yi 



Z/2 



-16Li2 — - 16Li2 — - 16Li2 - - 16Li2 - + IGLis - 



S2 



= -2sl + 2sl + 2sl - 8S2S3 + S2S3 



-321og'(yi) - 321og(?;3) \og{vi) - 321og(t;4) \og{vi) + 41og(vi) 
-8 log^(t;3) - 8 log^(t;4) - 4 log(t;2) - 4 log(t;3) - 16 log(t;3) log(t;4) 



-16Li2 



X2 



16Li2 



16Li2 



16Li2 



+ 16Li2 



2sl — 2s\ + 2^3 — 8siS2 + SiSs 

-32 \og\vi) - 32 log(t;3) log(^;i) - 32 log(t;4) log(t-i) + 4 log(t;i) 
-8 log^(t;3) - 8 log^(t^4) - 4 log(t;2) - 4 log(t;3) - 16 log(t;3) log(t;4) 



-16Li2 ( ^ 

Xi 



16Li2 ( ^ 

y2 



I6U, ( ^ 



16Li2 ( ^ 1 + 16Li2 



where the following notation was introduced: 



Xi = 




X2 = 


^{si-S2-Q^) + ^Vk, 


yi = 


^{si + Q^-S2) - ^Vk, 


y2 = 


^{S^ + Q^ -S2) + ^Vk, 


^1 = 


^{si + S2 + Q^)-^Vk, 


Z2 = 


^{si + S2 + Q^) + ^Vk, 


Vl = 




^^2 = 


7^(si-s2 + g') + -^v%, 

ZS2 ^S2 


^^3 = 




V4 = 




V5 = 




Vd = 


1-^. 



^2 



(25) 



(26) 



(27) 

(28) 
(29) 
(30) 
(31) 
(32) 
(33) 
(34) 
(35) 

(36) 
(37) 
(38) 
(39) 



We checked, that all infrared and ultraviolet divergences cancel as should be for vector interpolating 
currents. 

Now, let us proceed with the physical part of three-point sum rules. The connection to hadrons in 
the framework of QCD sum rules is obtained by matching the resulting QCD expressions of current 
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correlators with spectral representation, derived from a double dispersion relation at < 0: 

1 



— — — :^dsids2 + subtractions. 



(40) 



Assuming that the dispersion relation pn|l is well convergent, the physical spectral functions are 
generally saturated by the lowest lying hadronic states plus a continuum starting at some thresholds 
s^^ and s^2'- 



0{si 



e{s2 



(41) 



where 



9l 



2ml 



P^i{P2aP2f} +PlaPip) 



F,{Q')-F2m + 2{1 + 



2m2 



+P2pPlaPf.^ { 1 



2m^ 



F2iQ' 



1 + 



2m^ 



FsiQ' 



;P2aPlf3Pf,F3{Q^ 



2ml 



qMc^PiP - P2aP2p) [Fl(Q') + F2{Q^)] 



) [F^{Q') + F2{Q 



2ml 



+higher state contributions 



(42) 



The continuum of higher states is modeled by the perturbative absorptive part of H^ap, i-e. by 
Then, the expressions for the form factors Fi can be derived by equating the representations for 
three-point functions Hf^afB from ((7j) and ()40|) . It is reasonable to consider 3 sum rules: for structures 
Pitxgal3+P2fi9ci3, Pii,P2aPip + P2,iP2aPii3 and PipQaf, + P2c9 Ppi- This last step constitutes a formulation 
of QCD sum rules for our particular problem. 



3 Numerical analysis 



For numerical analysis we used Borel scheme of QCD sum rules. That is, to get rid of unknown 
subtraction terms in ((7j) we perform Borel transformation procedure in two variables si and 82- The 
Borel transform of three-point function n^a/3(si, S2, q') is defined as 



^l2n^a/3(Sl,S2,g^ 



lim 

n,m—> 00 



d 

dS2 



" „m+l 
■^1 



m\ 



Then Borel transformation 



^(pert|phys)^^2^^2^^2) 



of dZj) and (gO)) gives 

dsi I ds2 exp 
Jq 



(27r)2 



dsi 



si=mMf,S2=nM^ 



n/,a/3(Si,S2,g^) 



(43) 



S2 

Ml 



pt:T^'"'''\si.S2.q% (44) 



In what follows we put = M| = M^. If is chosen to be of order 1 GeV^, then the right hand 
side of (jl^ in the case of physical spectral density will be dominated by the lowest hadronic state 
contribution, while the higher state contribution will be suppressed. 
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Figure 3: dependence of electric form factor 



Now let us recall that our sum rules also receive power corrections proportional to QCD vacuum 
condensates. The evaluation of power corrections is simplified if performed directly for the Borel 
transformed expression of three-point correlation function and this is the reason why we delayed 
their discussion up to his moment. The quark condensate contribution is known already for a long 
time and is given by |2]: 

P/^iPlaPip + P2aP2f}) (13 - 18-^ + 18-^ + 2— 1 
+P,,P2f3Pla (10 + 36-^ - 36-^ + 4^ J + lQP^^Pl^3P2a 
+qf.{PlaPlf3 - P2aP2p) (5 - 18— + 18— + 2 — J 

HP2pgo., + Piagp,)M^ (64 + - 2^^ 

+b2P^g^pM^ + {pipgo.^.+P2.9p^)M^ (^-20 - 8-^^ |. (45) 

The correction due to gluon condensate was only partially computed in j2j, where subset of diagrams 
as well as part of Lorentz structures were considered. Taking everything into account we get^: 

= |^(oiG'p.G'^jo)| 

^See Appendix A for the details of calculation 
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4 2 / 4 

Pt,PipP2a-^ + PnPlaP2Pj^ + {pinP20P2a + P^nPl^Pla) ( 
-'2gij,fSP2a - '^PipQua " J^P2f39,j,a " J^Pla9fj.p + 5^2/35'/." + 5pia5'M/3 M46) 




Figure 4: dependence of F2 magnetic form factor 



Equating Borel transformed theoretical and physical parts of QCD sum rules we get 

,2, _ 9l fM'V 



= 3^ -T exp 



2m2 



dx' I dy's'six' + s's)A^^hicp[-x'] 



dx' I 
Jo Jo 

a., / ^ 



(47) 



5^ /M2^^ 



= ( — ) exp 



2ml 



30 



|(0) 



l(0) 



/"^ dx' / ^ dy' + 4)4°) - - - + 4)4"^ exp[-x'] 

^0 -'0 



+g / " / dy' [s[{y' + 4)4^) - 4(1/' - 4)4'^ - s',{x' + s',)A';> exp[-x'] 

^2 



2567r-^ 
'81M6 



a,(0|gg|0)^(4-4) 



3M4 ^ ' TT 



.10) , 



(48) 



= i ) exp 



87r^ yrrip 



2ml 



"0 



dx' I dy' 2s'^Af + 24 - 24^1^"^ - 8444^^ exp [-a;'] 
Jo 



1(0) 



1(0) 



l(0) 
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Figure 6: Borel mass dependence of p-meson magnetic form factor at = 1 GeV^ 

Here, for continuum subtraction we used so called "triangle" model. To verify the stability of our 
results with respect to the choice of continuum model we checked, that the usual "square" model 
gives similar predictions for pion electromagnetic form factor provided Sq ~ 1.5si is chosen^. In 

■^For more information about different continuum subtraction models see P] 
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what follows we set Sq = 2.2. This value coincides with one used in the previous analysis and 
is in agreement with the value of continuum threshold employed in the analysis of corresponding 
two-point QCD sum rules. For the rest of parameters entering our expressions for form factors we 
used the following values^: 

rup = 0.77 GeV, 



as 

71 

a^(0|gg|0)2 = 2.2 10-^ GeV^ 



{0\^GlJO) = 0.009 GeV^ 



In Fig. iniwe plotted the dependence of the p-meson magnetic form factor F2{Q'^) from the value of 
Borel parameter at = 1 GeV^. Similar dependence was also found for other form factors. To 
ensure stability of our estimates to the variation of Borel parameter we fix it value at = 2 GeV^, 
belonging to the " stability plateau" . 

The results for p-meson electric form factor including NLO corrections are shown in Fig. El (solid 
line is the sum of LO and NLO contributions, curve with long dashes denotes LO contribution and 
curve with short dashes stands for NLO contributions.). Similar results for magnetic F2{Q'^) and 
quadrupole F^IQ"^) form factors are shown in Figs. 0] and El correspondingly (again solid line is the 
sum of LO and NLO contributions, curve with long dashes denotes LO contribution and curve with 
short dashes stands for NLO contributions.). Unfortunately, our sum rules do not allow us determine 
magnetic moment of p-meson with precision better then already available predictions both from sum 
rules jSl HI and different quark models 013 CHI- At small our formula simply do not work - at 
this point we approach physical region of our three-point correlation function and OPE used in our 
analysis breaks down. Here one may only conclude, that the value of p-meson magnetic moment is 
close to 2 by extrapolating the ratio of F2{Q'^) / FiiQ"^) into the region of small momentum transfer 

Finally, let us following compare the behaviour of our form factors with radiative corrections 
taken into account in the limit Q"^ ^ oo with that, predicted by perturbative QCD (pQCD). The 
asymptotic of p-meson form factors is typically discussed in terms of states with given transverse 
and longitudinal polarizations of p-meson. In Breit system these two descriptions could be related 
with the help of following formula: 

FTTm = {p',eT\A,eT)/2E = F,{Q'), (50) 
FlAQ') = {p',eT\jf\p,eL)/2E = ^^iF^iQ') + F2iQ')), (51) 

Fll{q') = (p^eLbo^b,ei)/2i^^ = ^l(Q')-^^2(g^) + ^(l + ^)F3(g^), (52) 



where E is the meson energy E = yrrip + \Q'^. 

Quark counting and chirality conservation lead to the following asymptotic behavior of form 
factors: Fll{,Q'^) ~ 1/Q^, Flt^Q"^) ~ 1/Q^ and Ftt^Q"^) ~ 1/Q^- In terms of electric, magnetic and 
quadrupole form factors we have FiiQ"^) ~ F2{Q'^) ~ 1/Q^ and F^iQ"^) ~ It is easy to check 

that our form factors with radiative corrections included follow this behavior in asymptotic limit, 
while LO results contribute only as power corrections at large momentum transfers. 

*For numerical values of QCD condensates we took central values of estimates made in |S] 
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4 Conclusion 



We presented the results for p-meson electromagnetic form factors in the framework of three-point 
NLO QCD sum rules. The radiative corrections are sizeable (~ 30% in the case of Fi form factor 
and somewhat smaller for two other form factors) and should be taken into account when precision 
data on p-meson form factors become available. 

The work of V.B. was supported in part by Russian Foundation of Basic Research under grant 
01-02-16585, Russian Education Ministry grant E02-31-96, CRDF grant MO-011-0 and Dynasty 
foundation. The work of A.O. was supported by the National Science Foundation under grant PHY- 
0244853 and by the US Department of Energy under grant DE-FG02-96ER41005. 
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A 


Gluon condensate correction 



Here we present details on the evaluation of power corrections proportional to gluon condensate. 
This calculation could be relatively easy performed directly for the Borel transformed expression of 
three-point correlation function. Unfortunately, one of the methods (calculation in coordinate space), 
we will discuss below, does not allow us to subtract continuum contribution for gluon condensate 
correction. However, the form of the obtained expression leads us to the conclusion, that this 
contribution is simply absent in our final result. This conclusion is based on a fact, that typical 
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continuum contribution may show up only as incomplete F-functions. The latter are in fact present 
in contributions of each separate diagram, but they are canceling in the sum. 

The gluon condensate contribution to the three-point sum rules is given by diagrams with two 
external gluon vacuum fields, depicted in Fig. [7| For calculations we have used the Fock-Schwinger 
fixed point gauge [HI [T2j : 

x'A;{x) = 0, (53) 

where A^, a = {1, 2, . . . , 8} is the gluon field. The use of this gauge allows us express gauge potential 
An{x) in terms of field strength and its covariant derivatives at origin: 

= -l^uG^M - ^x,Xo^{D^G,,no) + (54) 

or in momentum representation: 

-\{2^Y {D^G.,{0)r gj^/ik) + ■■■ (55) 

So, basically, the calculation of gluon condensate correction is the ordinary calculation in back- 
ground of vacuum gluon fields in the form of ()54|1 or ()55|1 . Finally, vacuum averaging is performed 
according to rule: 

(0|G,.GJJ0) = ^5"^0|G^^G^^|0)(5,A. - (56) 

The calculation could be simplified if one uses the expression for quark propagator in background of 
gluon vacuum field: 

^i^^y) = -^^^"/5(o)w5 

\ 1 T If 

+ |^pj^^"^/3^"/5(0) - 1927^2 (^2)2 (^V - {xy f)Ga(3{0)Ga(3{0) 

+ operators of higher dimension, (57) 



where 



X — y, Gap — ^A^G^^, Go/3 — -^^aPuvG^y (5^ 



The corresponding momentum representation of quark propagator could be easily obtained through 
Fourier transformation. Next, for coordinates chosen as indicated in Fig. [7|from ()57p it follows that 
diagrams a) and d) are identically zero. 

Now there are two ways to proceed: first is to perform the whole calculation in momentum 
representation [2] and second way - do the same in coordinate space. Let's consider first calculation 
in momentum space. It is easy to find that the contribution of diagrams, where external gluon fields 
attached to different quark lines, is given by 



iV;3 (Pi,P2,gj - ^ 24 (27r)4 7 ^''ppfpi I 
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Figure 7: The gluon condensate contribution to three-point QCD sum rules. The directions of 
Pi, ki, k2 momenta are incoming, and that of p2 is outgoing. 
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' p2 p2 
-'1 -'2 

+ -p^Tr[7^Pi7"P27^fc + 27'^Pi7Vfc " (59) 

where Pi = pi + k and P2 = P2 + k. The integrals entering this expression could be conveniently 
evaluated using double spectral representation. For example, 

d'kk^^ = -.^-^ / d'k 



k^P^Pi 4 dpipdp2a J k^P^Pi 

- [ f do do ^"^^^^ - S2- Q^)PlM + -^1(^2 -Si- Q^)p2t, 

4:dpipdp2a J J ^ ^A3/2 (^si-pf){s2-p2) 

(60) 

and similar expressions hold for other two integrals. Finally, for the remained diagram c) the use of 
quark propagator in the background vacuum gluon field ()57|) allows present this contribution in the 
following compact form: 

7^ — ^ — 7^ — 7^ — S — n^a/3(Pi,P2,g), (61) 

_dpipdp2pdpiaP2a Opipdpip dp2adp2a \ 

where n^°2^(pi,p2, 9) is LO perturbative contribution to our three-point correlation function. Per- 
forming all differentiations and doing afterwards Borel transform according to 

s-(M=)(P^ = (^r^i-D-p^^-"'"' (62) 

we come to the final expression for gluon condensate contribution presented in the main body of the 
paper. 

Now, let us make a few comments about calculation of gluon condensate contribution within 
coordinate space representation^. The coordinate space amplitude corresponding to this contribution 
easily follows from an expression of quark propagator in the background gluon field ()57|). However, 
its Borel transformation is not that trivial. To do it, we first convert our result into momentum space 
with the help of the following formula [Tl] : 

_ / f^ip2X-ipiy a -la y ^ I ^) / grip2+^2p2_^T-3<?2 "^1 "^2 "^3 /ggN 

^4 y (^3, _ y-^2ly2m^2n 4'+"+™+l/!m!n! Jq {tiT2 + T2T3 + r3ri)'+'"+"-2 ^ 

The factors in numerator could be incorporated via: 

The Borel transformation of the resulting expression is performed with the help of the following 
formula: 

Bp2{M^)e~^P' = (5(1 - tM^) (65) 

The final expression for gluon condensate contribution obtained within this approach coincides with 
the result obtained in momentum representation and serves as a check of our result. 



^For more information see also appendix in 
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